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semiminor axis of an ellipse 

semimajor axis of an ellipse or semispan of a slab section 
surface pressure coefficient 

pressure coefficient behind oblique shock inclined 0° with respect to 
free-stream direction 

parameter defined by equation (3) 

length of cone 

free-stream Mach number 

virtual shock element (sketch (a)) 

element of actual shock 

surface pressure on cone plane of symmetry 
free-stream pressure 

pressure behind oblique shock inclined 0° to free-stream direction 
free-stream dynamic pressure 
corner radius of slab section 

length along body surface from nose apex (fig. 3) 
free-stream velocity 

average velocity normal to shock in shock layer 
axis tangent to oblique shock (sketch (a)) 
cone base reference length in X-direction 
axis normal to oblique shock (sketch (a)) 
angle of attack 

angle between shock and cone surface 
specific heat ratio of gas 

shock standoff distance for yawed cylinders 
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A s shock standoff distance on arbitrary body at location s 

AP pressure increment between shock and yawed cylinder surface in plane of 
symmetry 

Ax length increment of transformed body (sketch (a)) 
semiapex angle of cone in minor axis plane 
e z semiapex angle of cone in major axis plane 
0 inclination angle of shock to free-stream direction 
p 0 free-stream gas density 

Pq gas density behind oblique shock inclined 0° to free-stream direction 
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PREDICTION OF CENTERLINE SHOCK-LAYER THICKNESS AND PRESSURE 


DISTRIBUTION ON DELTA WING-BODY CONFIGURATIONS 
George E . Kaattari 
Ames Research Center 


SUMMARY 


Methods are presented to calculate both the shock inclination angle and 
the surface pressure coefficient in the vertical plane of symmetry of bodies 
at angle of attack. The methods are applicable over an angle-of-attack range 
from 0° to a maximum angle that depends on the body slenderness ratio; for 
very slender bodies, this maximum angle of attack approaches 90°. The 
methods apply to configurations of elliptical cross section and of rectangular 
cross section with rounded corners. 


INTRODUCTION 


There is much interest currently in pressures and shock- layer thicknesses 
generated on bodies at large angles of attack during atmosphere entry. Theo- 
retical solutions (refs. 1 and 2) have been developed for pressure distribu- 
tions and shock-layer thicknesses for symmetrical profiles at zero angle of 
attack. A recent work that uses numerical solutions of exact equations 
(ref. 3) has considered the case of nonsymmetric bodies at angle of attack. 
Approximate pressure distributions have been computed for more general (non- 
analytic) shapes such as wing-body combinations by tangent wedge and tangent 
circular cone methods and by Newtonian theory. 

The present method of estimating the shock locations and surface pressure 
coefficients in the vertical plane of symmetry is based on analysis of a sim- 
ple cone at angle of attack. The method may be extended, however, to the pre- 
diction of symmetry plane pressures and shock- layer thicknesses of configura- 
tions with arbitrarily varying (and expanding) cross sections. This extension 
to the method is accomplished by approximation of the body with a tandem 
series of cone frustums of appropriate cross section. In effect, this then is 
a tangent cone method but differs from the standard technique in that the tan- 
gent cone can be made to coincide closely with the entire windward local body 
periphery in the spanwise direction. This requirement is not possible with a 
circular cone that can be tangent only to a point on the body and thus cannot 
give a realistic shock-layer thickness. The present method, moreover, is more 
applicable to a larger angle-of-attack range than is the method of the tangent 
circular cone. 



ANALYSIS 


The analysis will be developed in two parts. The first part will discuss 
a semi empirical equation used to calculate the pressure coefficient in the 
vertical plane of symmetry on a conical body of arbitrary cross section at 
angle of attack. The equation is a function of the free-stream flow proper- 
ties, the shock inclination angle, 0, and the differential angle, 6, between 
the shock trace and the body surface in the vertical plane of symmetry. The 
second part presents a means of calculating the angle, 6, as a function of 
angle of attack, a, for conical bodies having two types of body cross section. 


Pressure Coefficient 

The pressure coefficient for a perfect gas behind an oblique shock 
angle 0 is 
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When the oblique shock is associated with a wedge (two-dimensional flow) , the 
angle between the wedge and shock is 0 = tan" 1 [ (p o /p 0 ) tan 0], where p o /p 0 
is the reciprocal shock density ratio across the oblique shock. The pressure 
coefficient on the wedge is identical to Cp (eq. (1)). If the shock is 

associated with an infinitely long (gy = 0) yawed cylinder, the shock trace is 
parallel to the cylinder and 0=0°. The pressure coefficient on the cylin- 
der is then due to the stagnation pressure corresponding to the normal Mach 
number Mq sin 0. The increase in pressure coefficient from the shock to the 
cylinder surface may be determined with the aid of shock tables. For normal 
Mach numbers, greater than 2.5, this increase is given to a good approximation 
(within 6 percent) by 
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It is assumed that the variation in the surface pressure coefficient 
between the limiting cases of a wedge and a yawed cylinder is a linear func- 
tion of the parameter k defined by 


tan 0 

(p O /p 0 ) tan 0 


(3) 


The "generalized" pressure coefficient is then 
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Equations (1), (2), and (3) when combined with equation (4) result in 

C = sin 2 0 - — — - sin 6 cos 0 tan 6 (5) 

p lY J Cy+DmJ 

When applied to elliptic cones at angle of attack, equation (5) was 
found to agree well with the results of a more rigorous analysis developed in 
reference 4 (see fig. 1). The centerline pressures given by the present 
method were also found to agree substantially with the values given by a semi- 
empirical technique presented in reference 5. 


Shock and Body Angle Relationships 

Equation (5) presumes a known relationship between the angles 0 and g. 
This relationship is not generally available from existing theories for coni- 
cal bodies of arbitrary cross section at angle of attack. This section, 
therefore, describes a method by which the shock and body angle relationships 
may be determined by simple calculations involving two-dimensional shock 
theory. 



Sketch (a) 

Sketch (a) depicts a conical body with its upper surface element in the 
vertical plane of symmetry defining an angle g with respect to the shock 
trace OX, which is inclined at the angle 0 with respect to the free-stream 
direction. If the horizontal component of velocity, V 0 cos 0, is neglected, 
the shock standoff distance A 0 from the base of the body at location x D to 
a "virtual" shock location, element m, is established by two-dimensional 
shock applied to the local body cross section in the normal flow field. When 
the horizontal component of velocity is considered, the standoff distance A 0 
is maintained with respect to an extension of the body surface element but is 
translated a distance Ax, so that the virtual shock element m appears at 
the real location n on the shock trace OX, this position being occupied 
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with respect to a transformed (elongated) body. Since the value, A Q , is 
established, 3 may be determined by the following procedure using the length 
of the transformed body x 0 +Ax. 


The Y-components of the stream velocity in the shock layer varies from 
the value immediately behind the shock, (p 0 /pe)V 0 sin 0, to the value 
approaching the body surface, V 0 cos 0 tan 6. The average velocity 

Vy. = (l/2)V 0 [(p 0 /pg)sin 9 + cos 0 tan 3]. The X-component of velocity is 

V Q cos 0. Therefore, a stream particle commencing from the shock at x Q is 
carried the horizontal distance Ax with velocity V Q cos 0 in the time 
interval required for it to approach the transformed body at Y-distance A Q 
with the average velocity Vy. Thus 



Ax V Q cos 0 


P0 


tan 0 + tan 3 


( 6 ) 


From sketch (a) 


A Q - x D tan 3 = Ax tan 3 


( 7 ) 


Substitution of equation (7) into equation (6) gives 


tan 3 ( tan 3 + — tan 6] + — /tan 3 - — tan 0) = 0 

\ p e I x o{ p 9 j 

in which x Q is related to the cone geometry by 

Z x o 

= — , b = Z tan c , a = Z tan e 

cos e cos 3 z y 
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(9) 


Substitution of the above value for x 0 into equation (8) gives the final 
result 


tan 3 [ tan 3 + — tan 0 

\ p e j 


A tan e cos e sec 3 
o z y 


tan 3 - — tan e ] = 0 

p e 


( 10 ) 

Equation (10) gives 3 as a function of 0; however, 3 as a function of 
angle of attack a is readily found through the relationship a = 0 - 3 - e 
(sketch (a)). ^ 
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The useful angle-of-attack range of the method is restricted in that when 
the centerline shock- trace incline approaches 0 = 90°, the angle 0 
approaches the approximate value (A Q /b)e z . The limiting angle of attack is 
then a * 90° - (A 0 /b)e z - e y . Thus, a reaches the highest value for slender 
cones (small e z and Zy ) , and for cones of a given slenderness, has higher 
values with increasing Mach number (decreasing A Q /b) . 

For small values of 0 (sec 0=1), equation (10) can be formed as a 
quadratic equation and an approximate solution readily obtained. The exact 
solution in chart form is presented in figure 2, wherein 0 is plotted as a 
function of (A 0 /b)tan e z cos z y for various parameter values (p 0 /p e )tan 0. 

The required value A Q /b is the shock solution for two-dimensional cylinders 
and is a function of the shock-density ratio and the body cross section. 
Solutions for various elliptic cross sections and slab sections with round 
corners are given in figure 3. These solutions are from an unpublished method 
that is a two-dimensional extension to the axisymmetric solutions of refer- 
ence 6. The discontinuities in the solutions are due to arbitrary step 
changes in the value of the gas specific-heat ratio, y, at Pe/p Q = 6 and 11. 
The effect of y may be estimated by interpolation if some "real gas" y can 
be specified. 


EXAMPLE APPLICATION 


The following numerical example illustrates the calculative procedure 
and use of the charts herein to determine the shock angle and pressure- 
coefficient solutions for a typical conical body at an angle of attack in air 
(y = 1.4) at Mach number 5. 

Given a/b = 1/3, z z = 15° (c y = 5.10°), and M Q = 5, determine the angle 
of attack and pressure coefficient when the shock is inclined (0 = 60°) to the 
free stream. The shock-density ratio, p 0 /pq, across the oblique shock with 
y = 1.4 is first determined: 

+ 2 2 q 2ii 

p e ~ (y+1) , 1Avi 2 . 2 . 2.4 2.4 x 25 x 0.750 U ““ A 

0 1 (y+1) M sm 0 


The reciprocal value, pq/p q =4.74. In figure 3(a), A Q /b = 0.692 at 
Pq/q 0 = 4.74 for an elliptic section with a/b = 1/3. 

Quantities required in figure 2 are now evaluated. The abscissa value 
(A 0 /b)tan e z cos e y = 0.692 x 0.268 x 0.996 = 0.185. The parameter 

(p 0 /pQ)tan 0 = 0.211 x 1.732 = 0.366 

The above values intersect on the curves of figure 2(a) at the ordinate value 
0 = 5.85°. The angle of attack 
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a = 6 - B - e =60°- 5.85° - 5.10° = 49.05° 

y 

The pressure coefficient (eq. (5)) is evaluated 

C = ill x 0.750 - - . 2 0 _ - 0.500 x 0.866 x 0.1025 = 1.297 
p 2.4 2.4 x 25 

Note that a specific angle-of-attack value cannot be selected a priori. 
Calculations with several values of 6 can be made quickly, however, and 
curves of Cp and 6 constructed as functions of angle of attack. The above 
procedure also applies to round-corner slab sections where r c takes the role 
of "a" (fig. 3(b)) . 


EXTENSION OF THE METHOD 


A delta wing-body configuration is depicted by solid outline in figure 3. 
The dashed lines represent outlines of a series of tandem cone frustums with 
which the configuration shape is approximated. The cone elements are con- 
structed tangent to the bottom surface at the vertical plane of symmetry and 
are also made tangent to the planform outline of the configuration as indi- 
cated. A typical cross section meeting these requirements is shown. Cross 
sections varying from a circle to ellipses of increasing ellipticity and 
several round-corner slab sections were used in the construction of figure 4. 
Thus, the actual configuration is approximated by a series of cone frustums 
of stepwise varying cross section, apex angle, and angle of attack. The 
accuracy of the approximation increases with the number of frustums used. 

Calculations of pressures and shock angles are then made for the individ- 
ual cones in the manner of the above example calculation. The shock- layer 
thickness is then determined by integration of the local shock angles (appro- 
priate to each cone) plotted as a function of body surface length: 

A = f S 6 ds = f S (dA/ds)ds (11) 

s 0 0 

Equation (11) may be integrated by either graphical or tabular methods. 


COMPARISON BETWEEN EXPERIMENTAL AND THEORETICAL RESULTS 


Cones 

Comparisons of experimental and predicted results of 8 as a function 
of a for elliptic cones are presented in figure 5. A wide range in cone 
geometry and Mach number (3 < Mq < 10) is represented. Agreement between 
experimental and predicted values is generally within 10 percent and is 
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particularly satisfactory at the higher Mach numbers. The data at Mach num- 
bers 2.94, 3.87, and 4.78 are previously unpublished results of tests in the 
Ames 1- by 3-Foot Wind Tunnel. 

A comparison between experimental and predicted centerline pressures for 
cones is presented in figure 6 in the form of a correlation curve. The ranges 
of angles of attack, Mach numbers, and configuration cross sections are noted 
on the figure. Experimental and predicted values correlate within the limits 
of about ± 2 percent. 


Wing-Body 

A comparison between experimental and predicted centerline shock- layer 
thickness for a delta wing-body configuration at M 0 = 7.4 is shown in figure 
7. Experimental and predicted values agree well. Note that the shock-layer 
thickness has the same value at a = 15° and 30°. 

Figure 8 presents a comparison between experimental and predicted center- 
line pressure distributions for the delta-wing and body configuration shown 
in figure 7. Agreement is satisfactory at angles of attack up to 30°. At 
a = 40° a discrepancy between experimental and predicted pressures occurs 
over a portion of the body length; however, the predicted values given by the 
present method are in better accord with experiment than are the indicated 
results of modified Newtonian theory. 


CONCLUDING REMARKS 


Simple methods were developed for prediction of inclination of a shock 
with respect to the windward trace on a conical surface in the vertical plane 
of symmetry and for estimation of pressure coefficient on this body trace. 

Predicted values were compared with experiments for air flows in the Mach 
number range 3 to 10 and the angle-of-attack range of 0° to about 70° for 
cones of a wide range in geometry. The methods were extended to the predic- 
tion of vertical symmetry plane shock-layer thicknesses and body surface pres- 
sures of configurations with arbitrarily changing cross sections with body 
length by the use of a tandem series of locally conical elements to approxi- 
mate the body shape. This, in effect, is a refined tangent cone procedure. 

The applicability of the methods has been demonstrated by comparison 
with test results from low enthalpy air flows. The suitability of the method 
for real gas flows is believed valid also since the real gas effect enters 
primarily through the density ratio across the shock. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, July 8, 1971 
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Figure 1.- Comparison of cone surface-pressure ratio given by the present 
method and that of reference 5, M n = 10, a/b = 2/3, e = 15.2°. 








0 .2 .4 .6 .8 1.0 


(A 0 /b) tan e z cos c y 

(b) (A 0 /b)tan cos £y = 0.2 to 1.0 
Figure 2.- Concluded. 
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Section A 


Figure 4.- Approximation of delta-wing-body by cone frustums. 
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Figure 7.- Comparison of measured and predicted centerline shock-layer thickness on delta-wing-body 
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